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Abstract.  In  this  paper,  we  present  the  results  from  Direct  Numerical 
Simulations  of  turbulent,  incompressible  flow  through  a square  duct,  with 
an  imposed  temperature  difference  between  two  opposite  walls,  while  the 
other  two  walls  are  assumed  perfectly  insulated.  The  mean  flow  is  sustained 
by  an  imposed,  mean  pressure  gradient.  The  most  interesting  feature,  char- 
acterizing this  geometry,  consists  in  the  presence  of  turbulence-sustained 
mean  secondary  motions  in  the  cross-flow  plane. 

In  this  study,  we  focus  on  weak  turbulence,  in  that  the  Reynolds  number, 
based  on  bulk  velocity  and  hydraulic  diameter,  is  about  4450.  Our  results 
indicate  that  secondary  motions  do  not  affect  dramatically  the  global  pa- 
rameters, like  friction  factor  and  Nusselt  number,  in  comparison  with  the 
plane-channel  flow.  This  issue  is  investigated  by  looking  at  the  distribution 
of  the  various  contributions  to  the  total  heat  flux,  with  particular  attention 
to  the  mean  convective  term,  which  does  not  appear  in  the  plane  channel 
flow. 

1.  INTRODUCTION 

In  this  work,  we  adopt  Direct  Numerical  Simulation  (DNS  hereafter)  as  a 
convenient  tool  for  studying  both  the  velocity  and  temperature  fields  in  the 
presence  of  solid  corners.  A square  duct  is  a suitable  geometry,  in  this  re- 
spect, due  to  its  high  degree  of  symmetry.  An  imposed  pressure  drop  drives 
the  flow,  while  the  temperature  field  is  generated  by  the  two  horizontal 
walls,  kept  at  constant,  different  temperatures  ±TW.  The  vertical  walls  are 
assumed  perfectly  insulated.  The  Reynolds  number,  based  on  hydraulic  di- 
ameter and  bulk  velocity,  is  approximately  4450,  while  the  Prandtl  number 
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is  assumed  0.71,  representative  of  air.  The  Reynolds  number  Rer,  based  on 
hydraulic  diameter  and  friction  velocity,  is  300,  as  in  [1]. 

Prom  the  present  simulations,  the  maximum  secondary  mean  velocity 
was  found  to  be  about  2%  of  the  centerline  streamwise  velocity.  The  results 
clearly  indicate  that,  even  though  both  the  heat  transfer  rate  and  friction 
factor  are  not  strongly  affected  by  the  presence  of  mean  secondary  motions, 
their  effect  on  the  distribution  of  various  quantities  within  the  flow  field  can 
not  be  neglected.  For  instance,  we  report  the  distributions  of  both  the  shear 
stress  and  the  heat  flux,  at  the  horizontal  walls.  Interesting  features,  clearly 
connected  with  the  presence  of  mean  secondary  motions,  characterize  these 
quantities. 


2.  PROBLEM  DEFINITION  AND  SCALES 


The  full,  time-dependent  incompressible  Navier  Stokes  and  energy  equa- 
tions are  solved  directly,  without  any  modelling  assumption  regarding  the 
turbulent  velocity  and  temperature  fields.  Simplifying  assumptions  are  ma- 
de, in  that  viscous  energy  dissipation  and  buoyancy  forces  are  neglected, 
and  the  fluid  properties  are  assumed  constant.  The  hydraulic  diameter, 
Dh  = 2 h,  the  friction  velocity,  uT  = (—APD/jApXx)1^2,  and  the  wall- 
temperature  Tw  are  used  in  order  to  define  dimensionless  quantities.  Do- 
main dimensions  are  indicated  in  figure  1,  which  shows  a sketch  of  the  duct. 
Two  simulations  have  been  considered,  and  they  differ  in  the  streamwise 
domain  length,  Ax,  which  was  47 xh  for  Case  A,  and  I2nh  for  case  B.  Grids 
with  200  x 127  x 127  and  600  x 127  x 127  computational  cells  were  used 
in  simulations  A and  B,  respectively,  so  that  the  resolution  was  identical  in 
both  cases.  The  grid  spacing  was  A:r+  = 9.42,  A y+  = A z+  = 0.45  -r  4.6. 
This  resolution  was  judged  to  be  adecquate  in  [1]. 

We  adopted  a second-order,  Finite- Volume  algorithm  in  order  to  solve 
the  flow  and  energy  equations  [2].  Decoupling  of  the  continuity  from  the 
momentum  equations  is  performed  by  a classical  second-order  projection 
scheme  [4].  The  computational  grid  is  uniform  in  the  homogeneous,  stream- 
wise  direction,  while  an  hyperbolic-tangent  distribution  is  adopted  in  both 
cross-stream  directions. 

All  data  reported  here  originate  from  case  B.  Statistics  were  acquired 
over  2 Large  Eddy  Turnover  Times  (LETOTS).  By  using  the  bulk- velocity 
reported  in  table  1,  this  time  interval  corresponds  to  about  1.6  flow-through 
times. 
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TABLE  1.  Global  statistics  for  Case  B, 
compared  with  results  for  the  channel 
flow,  and  with  a DNS  for  a square  duct. 


Case  B 

Channel  [6] 

Ref.  [1] 

Reb 

4453 

4560 

4410 

uh 

14.84 

15.2 

14.7 

f 

0.0363 

0.035 

0.037 

Nu 

4.89 

5.08 

3.  RESULTS  AND  DISCUSSION 

Mean  secondary  motions,  shown  in  figure  2(a),  result  in  four  couples  of 
counter-rotating  corner-vortices,  which  carry  high-momentum  fluid  from 
the  core-region  toward  the  corners,  and  low-momentum  fluid  away  from 
the  walls.  The  effect  of  the  mean  secondary  motions  on  the  streamwise 
velocity  distribution  is  evident,  in  that  isotach  lines  are  bent  toward  the 
corners.  Weaker  vortices,  located  on  both  sides  of  the  wall  bisectors,  were 
reported  in  [1],  but  have  not  been  observed  from  the  mean  data  of  the 
present  simulations.  A conditional  sampling  of  the  flow  field  showed  the 
presence  of  these  vortices  (figure  2b);  the  flow  field  was  acquired  only  when 
the  value  of  the  streamwise  vorticity,  averaged  in  the  marked  box,  was 
smaller  than  —3.0,  in  wall  units.  Since  these  vortices  are  highly  intermittent, 
a longer  averaging  time  would  be  required  to  capture  them  in  the  mean 
data.  At  higher  Reynolds  numbers  these  vortices  have  not  been  observed  [3], 
and  therefore  they  could  be  a low-Reynolds  feature. 

Mean  temperature  profiles  at  various  distances  from  the  left  wall  are 
presented  in  figure  3.  The  temperature  distribution  on  a wall-bisector  com- 
pares well  with  that  of  the  channel  flow,  at  the  same  ReT  [6].  The  tempera- 
ture profile  at  y+  = 35  from  the  left  wall  is  clearly  influenced  by  the  corner 
vortices,  resulting  in  a more  distorted  profile  than  at  the  wall-bisector.  The 
bending  of  the  temperature  profile  at  y+  = 3.3  indicates  the  presence  of 
significant  turbulent  activity  at  this  location.  This  is  quite  surprising,  since 
we  are  well  within  the  viscous  sublayer.  However,  the  side-walls  are  as- 
sumed perfectly  insulated,  thus  allowing  the  temperature  fluctuations  to 
extend  down  to  the  wall,  as  can  be  verified  in  figure  4(a).  Moreover,  the 
w'  fluctuations  assume  large  values  in  the  close  proximity  of  the  vertical 
walls,  before  vanishing  at  the  wall  (figure  4b).  Thus,  the  turbulent  mixing 
affects  the  mean  temperature  distribution,  even  in  the  viscous  sublayer  on 
the  vertical  walls. 

Values  of  the  bulk-Reynolds  number,  Re^ , the  mean  Nusselt  number, 
Nu  and  the  friction  factor,  /,  are  reported  in  table  1.  The  friction  factor 
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agrees,  to  better  than  1%,  with  empirical  correlations  available  for  this 
geometry  [5],  when  the  laminar  equivalent  diameter,  Di  = (l  + \/2j  h , is 
used.  The  Nusselt  number  is  defined  as  Nu  = (qw  2 h)  / [A;  (AT®)],  and  can 
be  evaluated  from  the  DNS  data  as 

1 / r+y 2 dT 

Nu  — - ( / RerPrwd  — — dy 

2 y— 1/2  dz 

where  the  brackets  indicate  avaraging  along  z,  and  all  quantities  are  dimen- 
sionless. The  Nusselt  number  from  present  simulations  is  only  3.8%  smaller 
than  the  Nusselt  number  evaluated  for  the  channel  flow  [6]. 

A possible  interpretation  of  this  fact  could  be  that  the  decrease  of  tur- 
bulent heat  transfer  ,near  the  vertical  walls,  is  partially  compensated  by 
the  mean  convective  heat  flux  in  the  core  region.  This  picture  is  confirmed 
by  figure  5,  which  shows  the  total,  the  turbulent  and  the  mean  convec- 
tive heat-flux  distributions.  The  gray  sfumature  corresponds  to  values  from 
— 11.9  (black)  to  11.9  (white).  The  total  heat  flux  ranges  from  —3.09  to 
11.90.  The  turbulent  heat  flux  assues  values  between  —4.5  x 10~5  and  6.24. 
For  comparison,  the  turbulent  heat  flux  in  the  channel  flow  [6]  assumes  an 
almost  constant  value  of  6.4  in  the  core  region,  when  scaled  by  the  fric- 
tion velocity  and  the  temperature  difference  between  the  walls.  The  mean 
convective  heat-flux  is  comprised  between  —6.05  and  9.95.  The  mean  con- 
vective heat-flux  is  defined  as  follows: 

<7conv  = 1/2  ReTPr  (W/uT)  (! T/Tw ) 

Figure  5 suggests  that  there  is  a corridor  through  which  the  heat  flows 
from  the  bottom  toward  the  top  wall.  Several  obstacles  tend  to  reduce 
this  flow.  First,  the  turbulent  heat-flux  is  reduced  near  the  vertical  walls, 
due  to  the  presence  of  the  velocity  boundary  layer.  Second,  the  corner 
vortices  create  mean  convective  heat  fluxes,  WT , which  pump  heat  in  the 
positive  2 direction  near  the  horizontal  walls  bisector,  but  they  also  create 
a negative,  mean  convective  heat  flux,  near  the  vertical  walls.  This  effect 
can  be  recognized  in  figure  5(c),  where  regions  of  negative  mean  convective 
heat-flux  appear  in  correspondence  of  the  corner  vortices. 

Figures  6 and  7 show  the  distributions  of  the  mean,  total  shear  stress 
and  heat  flux  at  the  bottom  wall,  respectively.  The  shear  stress  vanishes  at 
the  vertical  walls,  due  to  the  presence  of  stagnation  regions  at  the  corners. 
Three  local  maxima  are  evident  in  the  shear-stress  profile  at  the  wall.  The 
near-wall  peaks  are  located  in  the  region  below  the  corner-vortices,  and  are 
clearly  generated  by  the  intense  impingement  of  high-speed  fluid,  carried 
from  the  vortex,  against  the  wall.  The  local  minima  are  located  where  the 
corner  vortices  carry  low-speed  fluid  away  from  the  wall,  toward  the  core 
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region.  The  central  peak  in  the  wall  shear  stress  is  known  to  be  a low- 
Reynolds  number  effect  [1],  Clustering  and  widening  of  streamwise  velocity 
contours,  evident  in  figure  2,  allow  an  immediate  interpretation  of  the  shear- 
stress  behavior.  The  local  heat  flux,  figure  7,  shows  a behavior  similar  to  the 
shear  stress.  In  this  case,  however,  the  central  maximum  is  smaller  than  the 
side  ones.  The  lateral  peaks  in  the  wall  heat  flux  profile  are  located  below 
regions  of  high  mean  convective  heat  flux,  which  is  advected  by  the  corner 
vortices  and  carried  toward  the  bisector  of  the  horizontal  wall.  This  can  be 
appreciated  in  figure  5(c). 

The  w'9'  turbulent  heat  flux  decreases  near  the  vertical  walls,  due  to 
the  combined  effect  of  a reduction  in  the  intensity  of  both  w'  and  9',  as 
shown  in  figure  4.  Moreover,  there  is  also  a loss  of  correlation  between 
w'  and  9',  in  regions  close  to  the  vertical  walls,  as  reported  in  figure  8, 
which  shows  the  distribution  of  w'91  /w'9'.  In  figure  8 black  corresponds 
to  0 and  white  corresponds  to  1.  The  correlation  coefficient  is  particularly 
low  in  correspondence  of  the  corner  vortices  attached  to  the  vertical  walls 
{y/2h  « 0.4,  z/2h  & 0.35).  Just  below  this  area  there  is  a region  of  high 
correlation  coefficient,  attached  to  the  horizontal  wall.  It  corresponds  to  the 
high  heat  flux  region  evidenced  in  figures  5(a), (c).  The  loss  of  correlation 
of  w'  and  9'  close  to  the  vertical  walls  is  an  expected  feature,  considering 
the  different  boundary  conditions  imposed  on  velocity  and  temperature. 

Streamwise  power  spectra  for  both  the  temperature  and  the  streamwise 
velocity  fluctuations  are  shown  in  figure  9.  They  are  normalized  in  such 
a way  that  the  area  under  each  curve  equals  unity.  At  the  centerline  the 
spectral  distributions  of  temperature  and  velocity  fluctuations  are  very  close 
to  each  other,  in  the  whole  wavenumber  range.  In  the  corner  region  there  is 
still  substantial  agreement  between  9'  and  v!  spectra,  although  the  second 
tend  to  decay  slightly  faster.  As  one  might  expect,  in  the  corner  region 
both  viscosity  and  molecular  conductivity  are  more  effective  in  damping 
velocity  and  temperature  fluctuations,  and  both  spectra  decay  faster  than 
at  the  centerline,  at  least  in  the  intermediate  range  of  wavenumbers.  At  very 
large  wavenumbers  the  spectra  evaluated  at  both  locations  collapse  on  each 
other;  this  is  quite  surprising,  since  the  characteristics  of  both  the  velocity 
and  temperature  fields  are  fairly  different,  in  these  regions.  Aliasing  effects 
show  up  just  in  a very-high  wavenumber  range, which  accounts  for  negligible 
energy  content.  This  shows  that  the  adopted  streamwise  resolution  was 
adequate,  in  order  to  represent  the  main  features  of  both  the  flow  and 
temperature  fields. 

4.  CONCLUDING  REMARKS 

In  this  work  we  presented  results  obtained  from  DNS  of  turbulent  forced 
convection  in  a square  duct.  In  particular,  DNS  data  for  the  temperature 
field  were  not  available  before,  to  the  best  of  our  knowledge. 
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Streamwise  power  spectra  for  both  streamwise  velocity  and  temperature 
fluctuations  were  calculated,  at  various  cross-stream  locations.  Some  pecu- 
liar aspects  regarding  the  decay  of  the  highest  wave-number  components 
were  pointed  out. 

The  influence  of  mean  secondary  motions  on  both  the  temperature  and 
velocity  fields  was  discussed.  The  maximum  intensity  of  these  motions  was 
found  to  be  about  2%  of  the  maximum  mean  streamwise  velocity.  The 
comparison  of  the  duct  flow  with  the  plane-channel  flow  led  to  the  observa- 
tion that  global  parameters,  like  the  Fanning  friction-factor  and  the  Nusselt 
number,  are  not  strongly  affected  by  the  presence  of  the  side-walls,  while  the 
distribution  of  the  local  shear-stress  and  heat-flux  at  a wall  shows  charac- 
teristic patterns,  whose  origin  appears  to  be  related  to  the  mean  secondary 
motions. 

The  present  work  constitute  a base  for  future  investigations,  regarding 
unstably  stratified  internal  flows.  We  are  currently  running  other  simula- 
tions for  the  square-duct  flow,  in  order  to  investigate  the  combined  influence 
of  mean  secondary  motions  and  buoyancy  forces,  on  both  the  velocity  and 
temperature  fields. 
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Figure  5.  Heat-flux  in  the  2 direction,  (a):  total;  (b)  turbulent;  (c)  convective. 


